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Abstract 



The magnetized pure pion gradient (VV 9 ) phase in the deconfined Sakai-Sugimoto model is ex- 
plored at zero and finite temperature. We found that the temperature has very small effects on the 
phase. The thermodynamical properties of the phase shows that the excitations behave like a scalar 
solitonic free particles. By comparing the free energy of the pion gradient phase to the competing 
multiquark-pion gradient (MQ-VV 9 ) phase, it becomes apparent that the pure pion gradient is less 
thermodynamically preferred than the MQ-yy? phase. Bowever, in the parameter space where 
the baryonic chemical potential is smaller than the onset value of the multiquark, the dominating 
magnetized nuclear matter is the pion gradient phase. 
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I. INTRODUCTION 



Physics of dense nuclear matter is one of the most challenging area due to the lack of 
appropriate theoretical modeling. On one hand, the entities in the nuclear matter strongly 
couple to one another and therefore the perturbative treatment cannot be applied in a 
straightforward manner. On the other hand, the lattice approach to the Quantum Chromo- 
dynamics (QCD) can be applied to the situations of hot nuclear matter. The lattice results 
predict the deconfmement phase transition at temperature around 175 MeV for the dilute 
nuclear matter. However, this approach also faces difficulty in describing the nuclear matter 
with finite density due to the fermion sign problem. 

An alternative and complementary approach is the application of the holographic prin- 
ciple or the variations of the AdS/CFT correspondence l|-|3( in string theory to study the 
behaviour of nuclear matter. The Sakai-Sugimoto (SS) model 0, I^J is a holographic model 
which mimics the QCD at low energy most accurately. Starting with a type IIA string 
background with D4-branes as the source. Take the near-horizon limit and add the black 
hole horizon to generate Hawking-Page temperature to be identified with the temperature 
of the dual gauge matter. Since we need an approximately 4 dimensional QCD, one of the 5 
dimensional subspace is compactified into a circle whose radius is chosen so small that the 
Kaluza-Klein states are much heavier than the relevant energy scales and temperatures. 

The quarks and antiquarks are introduced as open-string excitations on the stack of 
U(Nf) flavour D8 and D8 branes located at fixed separation distance in the compactified 
coordinate. The boundary conditions of the sparticles in the circle are chosen to be anti- 
symmetric at the location of the flavour branes and the zeroth modes are thus eliminated. 
Consequently, the gauge theory at the flavour branes is a SUSY-broken 5 dimensional Yang- 
Mills theory with quarks and antiquarks in the fundamental representation. The effective 
theory has the same particle content as the QCD. Using the AdS/CFT correspondence, the 
bulk theory of this brane configuration is conjectured to be dual to the QCD-like gauge 
theory at the boundary. The striking feature of the SS model is that it provides a natural 
geometric realization of the chiral symmetry breaking. When the D8 and D8 merge at certain 
location in the radial coordinate, the quarks and antiquarks do not transform independently 
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under the chiral transformation and therefore the chiral symmetry is broken in the connected 
brane configuration. A chiral symmetric configuration occurs when the two flavour branes 
are parallel and the dual gauge matter will be in the chiral symmetric phase. 

Subsequent investigation reveals that the SS model accommodates the exoticpossibility 
that the chiral symmetry restoration and the deconfinement can occur separately 6J when the 
distance between the D8 and D8 branes in the compactified dimension is not too large. The 
deconfinement could occur at relatively low temperature while the chiral symmetry would 
be restored at larger temperature. Even though both the chiral symmetry breaking and the 
confinement are results of the strong coupling of the gauge theory, they are independent of 
one another as far as we know. It is thus possible that the real QCD also has distinctive 
chiral symmetry restoration and deconfinement. 

Chiral condensate of the QCD-like dense matter is explored in Ref. j3] using the Wess- 
Zumino-Witten induced anomalous term in the chiral perturbation theory and in Ref. jg] 
using the bottom-up AdS/QCD based on the confined SS model. When the magnetic field is 
applied, the condensate will respond by developing a gradient in the direction of the applied 
field. This gradient also carries the baryonic charge density proportional to the applied field 
and the gradient of the condensate. Holographic studies of the chiral condensate response to 
the magnetic field is investigated in Ref. [9| for the confined SS model. In Ref. 10], the pure 
pion gradient phase is explored and compared with the chiral symmetric quark-gluonplasma 
phase in the zero temperature approximation of the deconfined SS model. In Ref. Ill ], it is 
roughly compared with the mixed phase of the multiquark-pion gradient (MQ-y<^) using a 
zero-instanton limit of the multiquark configuration. The preliminary results suggest that 
the pure pion gradient phase might be thermodynamically less preferred than the MQ-yy? 
phase. In this article, we perform a thorough investigation into the pure pion gradient phase 
at finite temperature as well as its thermodynamical comparison to the MQ-y<^ in order to 
obtain a more definitive quantitative result. It is found that the pure pion gradient phase is 
insensitive to the change of temperature in the range T = — 0.16. It is also shown that the 
pure pion gradient phase is generically less preferred than the MQ-y<^ phase except when 
the baryon chemical potential is smaller than the onset value of the multiquarks. In that 
region of the phase diagram, the dominating phase is the pure pion gradient. 
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The article is organized as the following. In Section II, we setup the holographic model 
of the magnetized chirally broken nuclear phase without an instanton. A zero temperature 
solution is obtained and relevant dual physical quantities as well as their relationships are 
discussed. Thermo dynamical properties of the pure pion gradient phase at finite temperature 
and the comparison with the multiquark phase are discussed in Section III. Section IV 
concludes the article. 



II. HOLOGRAPHIC SETUP OF THE MAGNETIZED CHIRALLY BROKEN 
PHASE 

In the non-antipodal SS model, a stack of N c D4-branes generates a curved 10 dimensional 
spacetime in type IIA string theory. The near-horizon limit of this background is then taken 
and the black hole horizon is added by introducing the factor f(u) [s,[l2] into the background. 
The x 4 direction is compactified with certain radius to obtain an effective (1 + 3) dimensional 
subspace in the low energy limit. The resulting spacetime of the Sakai-Sugimoto model is in 
the form 

ds i = (^-j 12 (f(u)dt 2 + 5 l3 dx l dxi + dx A 2 ) + (u 2 dn 2 4 + j^j 

e = g s I £J , R 3 m = *9sNJ 3 s , 

where f(u) = 1 — Uj./u 3 , ut = 167r 2 i?£, 4 T 2 /9. T is the Hawking-Page temperature of the 
black hole which is identified with the temperature of the dual gauge matter at the boundary. 
Rr>4 is the curvature of the background which is generically different from the compactified 
radius R of the x 4 coordinate. is the dilaton field, a function of u in this background. 

We then introduce stacks of Nf D8 and D8 flavour branes with separation L on the circle 
of compactified x 4 at the boundary u — > oo. Open string excitations with one end on these 
branes behave like chiral "quarks" and "antiquarks" in the fundamental representation of the 
U(Nf). In the brane configuration where D8 and D8 are parallel, open-string excitations on 
each stack of branes transform independently under the chiral transformation and thus we 
have a chiral symmetric background. The dual gauge matter will be in the chiral symmetric 



phase. On the other hand, in the connecting brane configuration, chiral symmetry is broken 
at the tip and the corresponding gauge matter will be in the chirally broken phase 6]. 

To add the baryonic density to the boundary gauge matter, the non-normalizable mode 
of the (Iq component of the U(l) C U(Nf) field is turned on. The baryon chemical potential 
fi of the corresponding gauge matter is identified with the non-normalizable mode of the 
DBI gauge field at the boundary by 13] 

[l = clq(u — > oo). (1) 

To turn on the magnetic field, another component a\ is used as the vector potential 
generating the magnetic field. The direction of the magnetic field is chosen so that the 
vector potential is 

a\ = Bx 2 . (2) 

The Chern-Simon action in the background couples these two components to the third 
component a A of the U(l), generating the response to the external magnetic field. The 
response appears as the gradient of the chiral condensate along the direction of B at the 
boundary which is defined to be a A {u — > oo) = yyj. Here and henceforth, we will call \/ip a 
pion gradient. 

The DBI and the Chern-Simon actions are then given by 



Sds = du u^ 2 Jl + + f(u)(a' A ) 2 - «) 2 + f(u)u^x' 2 , (3) 

o POO 

Scs = -zM / du {d 2 ala v af - 9 2 a 3 y a V), (4) 

^ Juc 

where M = N c R 2 Di / '(67r 2 (27ra') 3 ) defines the brane tension. The factor 3/2 in the Chern- 
Simon action comes from addition of surface term in order to maintain the gauge invariance 
of the total action in the situation when the gauge transformation does not vanish at the 
boundary (see Ref. {9] for details). 

Consequently, the equations of motion with respect to each gauge field , are 



Vu 5 + B 2 u 2 f(u)a[ A . 3 



j A --Bfi + 3Ba%, (5) 



y/l + f(u)( a[ A y - K ) 2 + 2 

Vu 5 + B 2 u 2 a' y , 3 . o . 

= = d Ba A (00) + 3Ba? . (6) 
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The corresponding density and current density, d, Ja, at the boundary(w — > oo) are defined 

as 

5S eo -i 



= (d,f A )- 



They are related to the components of the U(l) gauge field by 

Vu 5 + B 2 u 2 < 



J-4 



Vi + /H«) 2 - KT) 2 + f{u)u^ 



Vu 5 + 5 2 u 2 /(«)aj 



(7) 
(8) 

(9) 
(10) 



Vl + /(«)(a^) 2 -«) 2 + /(«)«3x? 
For the phase of pure pion gradient where chiral symmetry is broken, the axial current ]a 
is set to zero and the density d = ^Bxjip is the definition adapted from the Wess-Zumino- 
Witten action of the boundary gauge theory 

The constant of motion with respect to x^{u) for the pure pion gradient phase yields 

D(u) 2 ) 



{<{u)Y 



i 



u 3 f(u) 



u 3 [f(u)(C(u) 



9B 2 (a%-!± 



F 2 



where 



D(u c ) 2 ) 



9B 2 [aV(u c ) 



2 ) ^4(^0 J 



y/l + f( Uc )ul X' 2 {U C ) 



ul'\ f(u c )C(u c 



9B 2 [a%{u c 



(12) 



(13) 



where C(u) = u 5 + B 2 u 2 , D{u) = d — 3B\/(p/2 + 3Baf(u). u c is the position where the D8 
and D8 branes connect. Since there is no instanton in this case, the branes connect smoothly 
at u c . We also have D{u c ) = from a^(w c ) = 0, and x' A {u c ) = 00. 

Since the DBI action, Eqn. ([3]), is divergent from the limit u — > 00, we would need the 
action of the magnetized vacuum for the regularization. For the magnetized vacuum, we can 
let the non-normalizable modes, a^,af = and d,j A = 0. The vacuum action then takes 
the following form 



S [magnetized vacuum] 



yJC(u)(l + f(u)u*X%] 



du. 



where 



x' 4 {u)\ vac = 1 =. (14) 



\jn u ) u {f(u Q )uic( U0 ) 



Again, the position of the tip of the brane configuration is denoted 



by u . The temperature 



and field dependence of the position uq are given in Fig. 1 of Ref. [14j. It saturates approx- 
imately at 1.23 for all temperatures at high magnetic field. The action of the vacuum will 
be used to regulate the infinity of the DBI action from the limit u — > oo when we calculate 
the free energy of the dual gauge matter in the subsequent section. 



A. Zero Temperature case f(u) = l 

We can numerically solve the equations of motion, Eqn. (EJ) , (J6j) by using the shooting 
algorithm. However, it is illustrative to consider first the limiting case of zero temperature 
where f(u) = 1 and the equations of motion are simplified so much that they yield exact 
analytic solutions. Later on we will actually find from the numerical solutions that most 
physical properties of the pion gradient phase are insensitive to the change of temperature. 

In the zero temperature case, starting from the equations of motion, Eqn. ©,([6]), rewrit- 
ten as 

- £) = f(u)afa[ A (15) 



2 

Qf>2 A( V _ l±\ 

f(uKa? = s Ql lQ ° a * J (l + « + K 42 " O- (16) 



From Eqn. (TIB]) , for / = 1 we can solve to obtain 



i4 - \f - (^) 2 = af - (W) 2 . (17) 



Using Eqn. ( ITbl) . direct integration leads to 



V W) 2 cosh/(w) (It 



4 = a/(|) 2 -(V^) 2 sinh/( M ), (19) 
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where 



I(u) 
g(u,u c ,B) 



du\ g(u,u c ,B) 



C(u) 
9B 2 



V\2 



i + u 3 xf 



1 + 



n 3 (CH-9i? 2 ((f) 2 -(v^) 2 )) 
95 2 ((^ 2 



-i -1 



(20) 
(21) 
(22) 



_n3(CK)-95 2 ((f) 2 -(v^) 2 )) 
by using the boundary conditions af(u c ) = and Eqn. (fTTj) at it c . Additionally, there are 
two constraints which need to be satisfied, 

fjt 1 



where 

The pion gradient is thus 



cosh/^, 

2 x ^ 

1 = 2/ rfu 2:4(11), 

Ac 

I(it — > 00). In the zero temperature case x 4 is given by 
. n3(C(n)-9fi 2 ((f) 2 -( V ^) 2 )) 



(23) 
(24) 



u ' 



M 3(CK)-9 J B 2 ((f) 2 -( W ) 2 )) 



-1/2 



— tanh J n 



(25) 



(26) 



In order to obtain the solutions, we numerically solve for u c from the constraints Eqn. ( 123]) 
and Eqn. ( 124)) simultaneously by fixing two parameters among (B,d,^/cp, fi). The solutions 
always have | > VV 9 as a reality condition. 



As is found in Ref. [10J for the pure pion gradient and Ref. [ill, [l_5| for the model with 
instantons, there are 2 possible brane configurations satisfying the scale fixing condition 
Lq = 1, one with small and one with large u c . The brane configuration with small u c has 
longer stretch in the w-direction and therefore has higher energy than the configuration with 
large u c . The excess energy makes this configuration less preferred thermodynamically. For 
the pure pion gradient phase, there is also the small-w c configuration (for sufficiently large \x 
and small B) which we found to be less preferred thermodynamically even than the vacuum. 
Therefore this configuration will not be considered in this article. 

The solutions can be explored by slicing through the plane in the parameter space at 
fixed magnetic field (B), fixed density (d), and fixed chemical potential (/j,) respectively. At 
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fixed B, the solutions are shown in Fig. [TJfor the position u c and the chemical potential as 
a function of the density. For small B, the position u c has certain variation with respect to 
the density. As B increases, u c saturates to an almost constant curve with a slight density 
dependence. The chemical potential at fixed B(> 0.1) is found to be an exact linear function 
of the density. The slope of the linear function is inversely proportional to B. The relation 
can be summarized into the following simple form 

H = ~ for B > 0.1. (27) 

This implies from d = 3B\/ip/2 that fi = 2\/(p for B > 0.1. The linear relation between /i and 
d can be interpreted as the absence of self-interaction among the pion gradient excitations. 
Each pion gradient excitation seems to behave as free entity for B > 0.1. 



u c fl 




01234 01234 



FIG. 1: The position u c (a) and the chemical potential (b) as a function of the density at fixed 
magnetic field B at zero temperature. 

For fixed d, the position u c and the chemical potential are shown as functions of B in 
Fig. [2|3l for d = 1.0. Solutions exist for the entire range of £>, down to u c {B = 0) = 0. 
Consistent with what previously found, \i is found to be inversely proportional to B as is 
shown in Fig. El 

For fixed /x, Fig. H] shows interesting transition between 2 regions of the parameter space. 
In Fig. El the relation between d and B is shown to be approximately quadratic for B < 0.2 
and linear for B > 0.2. From d = 3B\/(f/2, this implies that VV 9 is a linear function of B for 
B < 0.2 and a constant function for B > 0.2. Fig. S(b) confirms the behaviour. Since the 
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u c,0 




0.5 1.0 1.5 2.0 2.5 3.0 



B 



FIG. 2: Position u c as a function of B at a fixed density d = 1.0 at zero temperature, the position 
uq of the magnetized vacuum is shown for comparison. 




slope = -0.999 
-4 -3 -2 



log(fl) 



FIG. 3: Chemical potential as a function of B at d = 1.0, T = in the logarithmic scale. 

saturation at large B occurs around VV 9 = A*/2, the slope of the linear region, B < 0.2, is 
therefore proportional to \i. Consequently, for small B, xjip ~ fiB. The behaviour at small 
B is similar to the behaviour found in Ref. jl, 9] for the confined phase. The result in the 
deconfined phase of the SS model at zero temperature was first obtained in Ref. 10]. 

It should be noted that the nonlinear effects of the DBI action become apparent for 
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B 



0.4 
0.0 0.5 



FIG. 4: The density, the pion gradient, and the position u c as a function of B at fixed fx = 0.2, T = 0. 



> 0.2 where VV 9 — / i /2- As B increases, the pion gradient does not change but its baryonic 
density increases linearly with the field. This y '^-saturation is a new effect observed only 
in the theory with DBI gauge interaction. 



III. THERMODYNAMIC AL PROPERTIES OF THE PURE PION GRADIENT 
PHASE 

In the pure pion gradient phase, since x' A — > oo at u c , the integrand of the action diverges 
at u c in addition to the limit u — > oo. This also occurs with the magnetized vacuum where 
x' A is divergent at uq. However, the limit which makes the integral and consequently the 
action divergent comes only from u — > oo (the divergences at wo, c are weaker than a simple 
pole and thus finite over integration). We can therefore regulate the action by subtracting 
the total action with the action of the magnetized vacuum in the usual manner. 



II 



logW 



// — 9 — > 




/i 

—4 




-6 




-3 -2 -1 


1 



log(5) 



FIG. 5: The density as a function of B at fi = 0.2, T = in the logarithmic scale. 

The regulated free energy is thus given by 

jF E (d, B) = B) + /id, (28) 

where Q(^,B) = S[ao(u), ax(u)](e.o.m.) — S [magnetized vacuum] = J-{n, B). 

We can calculate the total action satisfying the equation of motion 
S[a (u), ai(u)](e.o.m.) = Sc 8 + Scs to be 



,5' 



D8 



N / du C(u) 



f(u)(l + f(u)u*x?) 



\ f(u)(C(u) + D(uy)-9B2(aX-^ 



2 ■ 



s, 



cs 



o poo 

-M-B / d% 



3J3&Q (&Q 2 ; 



f(u)D(u)af)Jj^ 



f(u) 1 " X 4 



f(u)(C(u) + D(uy)-9B* 



(29) 



(30) 



For zero temperature the total action reduces to 



r>oo 

J\f I du\ 



g(u,u c ,B) 



C(u)-9B 2 ((^ 2 



C{u) 



95 V v 



2 v 2 



(w) 2 ) 

We can compute this action by substituting Eqn. ( Tl8|) into the expression. The free energy 
at fixed chemical potential J-"(/i, B) of the pure pion gradient phase at zero temperature 
is shown in Fig. [6j Once d, /i > 0, the free energy becomes smaller than the free energy 
of the magnetized vacuum (being negative) and thus thermodynamically preferred. The 
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magnetization at fixed chemical potential B) = — d ~ F ^ I} B ^ therefore increases from zero 

and becomes constant M(/x = 0.2, B) ~ 0.0152 at large field (B > 0.2) as we can see from the 
slope of Fig. El On the other hand, the magnetization at fixed d — 1.0, M(d, B) - 



dF E {d,B) 



dB 



of the pure pion gradient phase is a rapidly decreasing function of B as is shown in Fig. [7J 
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FIG. 6: The free energy as a function of B at fixed [i = 0.2, T = 0. 
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FIG. 7: The free energy as a function of the density of the pure pion gradient phase compared to 
the multiquark-yi/? phase at fixed d = 1.0, T = 0. 



The pressure of the pure pion gradient as a function of the density can be calculated 
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using Eqn. ([27]) and d = % (see Ref. (l6|) 



dfi 

P(d,B) = n(d,B)d- [ fi(d',B) d(d'), (31) 

Jo 

= \KB)d\ (32) 

where k(B) = A/3B for B > 0.1. The quadratic dependence of the pressure on the density 
without higher order term reveals that the pion gradient excitations behave like free particles 
without either repulsive or attractive interaction among themselves. 

For the parameter space in the region dCl,.B<l such as the regions shown in Fig. H] 
13 since yy? ~ pB,d = 3B\/ip/2, we have d = apB 2 for some constant a. In this case, 
the linear relations between p and d is still valid and the equation of state is again given 
by Eqn. ( 132]) with k(B) = 1/aB 2 (for p = 0.2, a ~ 4.634). This behaviour is similar to 
what found in Ref. [7[ using the Wess-Zumino-Witten term in the boundary theory and in 
Ref. jsl [hJ] for the confined and deconfined SS model at zero temperature. 

The energy density can be calculated straightforwardly 

P = I H(V, B ) dr], (33) 
= \KB)d\ (34) 

where k(B) = l/aB 2 ,4/3B for small and large B respectively. The results are remarkably 
similar to the results from the bottom- up AdS/QCD model considered in Ref. jg]. The 
equation of state then becomes simply P = p representing free gas of the solitonic excitations 
of the the pion gradient. The adiabatic index, T, and the sound speed, c s , are then calculated 
to be 

^dP _ p 2 _ 

p dp ~ p Cs ~ ' [6b) 

c s = 1, (36) 

the typical behaviour of the free gas. 

For nonzero temperature, the full equations of motion, Eqn. (E]),® can be solved nu- 
merically by double shooting algorithm aiming for two conditions to be satisfied at once: 
af(u c ) = 0,L = 1 while fixing B,p and d (and consequently VV 9 )- It is found that the 
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temperature dependence of every physical quantity of the pure pion gradient is very weak. 
Figure [8] shows the difference of the density, the pion gradient, and the position u c at fixed 
fi = 0.2 between T = 0.14 and T = 0. Observe that the difference in the temperature 
dependence of the density and the pion gradient are the most distinctive in the transition 
region when the magnetic field changes from small to large values. 
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FIG. 8: The difference between the density, the pion gradient, and the position u c of the yy? phase 
at T = 0.14 and T = as a function of B for u = 0.2. 



A. Comparison to the Multiquark-y^ phase 



We would like to consider whether the pure pion gradient phase is thermodynamically 
preferred than the other possible nuclear phases exist in the same parameter space. In the 
deconfined SS model in the presence of the magnetic field, there are generically 3 possible 
phases in addition to the vacuum; the chiral symmetric QGP phase, the chirally broken phase 
of multiquark-y</? (MQ-yyj), an d the pure pion gradient (yy?) phase. The multiquark 



nuclear phase has been studied in Ref. Ill, ll4J, [17j and found to be the most preferred 
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phase for the dense deconfined nuclear matter under moderate external magnetic fields. The 
multiquark phase actually has certain mixture of the pion gradient as the source for the 
baryon density. This is inevitable since the response of the nuclear matter to the external 
magnetic field is in the form of the spatial variation of the chiral condensate in the direction 
of the applied field which we call the pion gradient. 

However, the ratio of the pion gradient population with respect to the multiquark de- 
creases as d grows jl4j . It is thus suggestive that the multiquark phase is likely to be more 
thermodynamically preferred than the pion gradient phase. In this subsection we directly 
compare the two phases at zero temperature using the free energy at fixed density d — 1.0. 
The MQ-\/(p phase imposes the boundary conditions (see Ref. [n], Q for details); 

3a = 0, af(u c ) = 0, a%{u c ) = -u c ^/ f(u c ) + n s (u c - u T ), 

where n s is the number of colour strings (hanging from the baryon vertex down to the 
horizon) in fractions of 1/N C . 

The result is shown in Fig. [7J Clearly, the multiquark-pion gradient (MQ-yy?) phase 
is more preferred than the pure pion gradient phase. Similar behaviours are confirmed for 
small d ~ 0.1 and large d ^> 1. Especially at large densities, since the baryon chemical 
potential of the MQ- \/ip phase increases slower than a linear function ll| whilst it is linear 
for the x/(p phase, the dominant term [id in the free energy for the MQ-y<£> phase becomes 
much smaller and thus more stable thermodynamically. 

However, there is a region of parameter space where the yy? phase is dominant. When 



the baryon chemical potential < fi ODSCt = [±{d = 0) = |n 0A //(u ) + n s (u — u T ) of the 
multiquark, the multiquarks cease to exist and the pion gradient which can be constructed 
at arbitrarily small \x (since \i ~ d) will be dominating. The corresponding transition line in 
the (fi,T) diagram for B = is shown in Fig. 8 of Ref. 17j. For B > 0, dependence of Uq 
on B affects the transition line accordingly as shown in Fig. [9j The dotted line represents 
schematic transition to the chiral symmetric quark-gluon plasma (xS-QGP) phase. The 
chiral symmetry restoration between the magnetized vacuum and the x^-QGP has been 



studied in Ref. 



The transition between the pure pion gradient phase and the x^- 



QGP has been explored in Ref. 10] with / = 1 approximation for the pure pion gradient. 
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FIG. 9: The onset chemical potential of the multiquark- xjip phase as a function of T, B (for 
B —7- oo, tio = 1-23 is used). These lines can be served as the transition lines between the \/(p phase 
on the left and multiquark-vv 3 phase (n s = 0) on the right. The dotted line represents schematic 
transition to the chiral symmetric QGP phase. 

Since we found that the VV 7 phase is insensitive to the change of temperature, the results in 
Ref. 10] should be justified to be a good approximated phase diagram. The chiral symmetry 
restoration between the MQ-yv 9 phase and the X'S'-QGP phase has been investigated in 
Ref. Q. 



IV. CONCLUSIONS AND DISCUSSIONS 



The behaviour of the chirally broken pure pion gradient phase in the deconfined SS model 
is studied at zero temperature and subsequently finite temperature. The magnetic response 
of the chirally broken phase is linear \/ip ~ fiB for small field and saturates to constant 
value VV 9 ~ /x/2 for large field. 

Relationship between fj, and d is also linear \x = k(B)d where k(B) ~ 1/B 2 , 1/B for small 
and large B respectively. This implies that the excitations of the pion gradient behave like 
a free gas with no interaction among each other. The equation of state is thus simply P = p 
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with the sound speed equal to the speed of light. The free energies at fixed fi and d are 
obtained numerically. Magnetization at fixed \i increases with B for small field and drops to 
constant value for large field. Magnetization at fixed d is a decreasing function with respect 
to the magnetic field. 

Using the free energy at fixed density, we show that the pure VV 9 phase is less preferred 
thermodynamically than the MQ-y</> phase at zero temperature. The configuration of the 
pure pion gradient phase is found to be insensitive to the change of temperature, the differ- 
ence of the free energy at fixed /i for T = and T = 0.14 is minimal, only about < 2x 10~ 4 . 
On the other hand, the free energy of the MQ- y<^ phase is a decreasing function in the tem- 
perature Jjj]. Therefore, we can conclude that the pure pion gradient phase is generically 
less preferred than the MQ-\/ip for general situation. 

However, there is an exception for small chemical potential, \x ~ 0.175 — 0.41. When 
A* < AWsct of the multiquarks, the multiquarks simply cannot exist while the pion gradient 
can be induced at arbitrarily small \i. Therefore, in this region of the parameter space, the 
pure pion gradient phase is dominating over any other phases. The transition lines are given 
by fi = /ionset in the (/i, T) plane. The interior of certain classes of the dense astrophysical 
objects such as the maguetars Q would have the correspond* g regions where the chemical 
potential (and the density) and temperature fall into this range. In those regions, the 
dominating nuclear phase which governs physics of the stars would be the pure pion gradient 
phase. 
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